Oscillating Ponomarenko dynamo in the highly conducting limit 
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This paper considers dynamo action in smooth helical flows in cylindrical geometry, oth- 
erwise known as Ponomarenko dynamos, with periodic time dependence. An asymptotic 
framework is developed that gives growth rates and frequencies in the highly conducting 
limit of large magnetic Reynolds number, when modes tend to be localized on resonant 
stream surfaces. This theory is validated by means of numerical simulations. 



I. INTRODUCTION 

A well-known kinematic dynamo model goes back to the work of Ponomarenko who found 
that magnetic modes could be amplified in a flow field in cylindrical geometry (depending only 
on distance from the axis), which generally possesses helical streamlines. In recent studies this 
Ponomarenko dynamo has been investigated when the helical flow is modulated in time 0,0], with 
a focus on the dynamo threshold. The aim here is to quantify the effect of simple time-periodic 
fluctuations on the mean flow, and the effect of these on the threshold for magnetic growth. The 
main conclusion is that the dynamo threshold is larger than the one obtained without fluctuations, 
suggesting that large scale fluctuations are not desirable when optimizing a dynamo experiment. 
In dynamo experiments, such large scale fluctuations have been avoided simply by adding inner 
walls [3) 0] or a flow-stabilizing ring 0]. A further experiment in preparation [?J is based on a 
single helical flow, again avoiding large scale fluctuations. 

The above simulations 0, 0| were purely numerical and in order to give some theoretical backing 
to the results it is necessary to use an asymptotic limit where approximate results can be obtained. 
The steady Ponomarenko dynamo [l| has been studied for Rm S> 1 in the kinematic case 0, 0, Q3] 
and equilibrated solutions have been found taking into account the nonlinear feedback on the flow 



111 ]. In both cases the underlying flow is steady, and our aim here is to extend the kinematic 
results to the case of oscillatory flow fields. 

We therefore adopt the limit of large magnetic Reynolds number (Rm 3> 1), generally much 
above the threshold. Although our aim is to derive asymptotic results which are valid for this 
general class of flows, with an eye to experimental dynamos and numerical simulations, we note 
that Ponomarenko dynamos may also occur in bipolar jet-like outflows commonly observed in 
protostellar systems, in which magnetic field probably plays an important role. Though it is 



argued that such magnetic fields are produced inside the protostellar disk [12|] we cannot exclude 
the existence of a Ponomarenko type dynamo in such a helical jet in which strong time fluctuations 
may also occur. 

Here we assume a time-periodic flow depending only on radius, limit our investigation to the 
kinematic approximation, and study the asymptotic limit of large Rm (section [TT]) . Our analysis 
will be compared to direct numerical simulations for three cases (section llllj) : stationary flow and 
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oscillatory flow with zero mean flow (ZM), or non-zero mean flow (NZM). Finally in section HVl we 
generalize our analysis to the case of time- varying resonant radius. 

II. MODEL AND ASYMPTOTIC APPROXIMATION 

The time evolution of the magnetic field is given by the dimensionless induction equation 

e— = V x (U x B) + eV 2 B, (f ) 

with e = Rm" 1 and where we have adopted a diffusion time scale for our time variable (in contrast 
to 0], but in accord with 0]). We consider a time-dependent helical flow expressed in cylindrical 
coordinates (r, 9,z) by 

U(r,t) = (0,rn(r),V(r))F(t) for r < 1, U = for r > 1, (2) 

where Q and V are smooth functions of r and F is a given function of time. In the stationary 
case F = 1, and for the non-stationary case we will consider two functional forms for the time- 
dependence, 

F(t)=cosujt (ZM), F(t) = I + p cos cot (NZM), (3) 

the 'zero-mean' and 'non-zero mean' flows, respectively. 

For the linear, kinematic dynamo problem we may consider a magnetic field of the form 

B (r, t) = e ^ m9+kz+ ^ b(r, t) , (4) 

where m and k are the azimuthal and vertical wave numbers of the field and (ft(t) is a phase, put 
in for convenience, that we will choose shortly. The solenoidality of the field, V • B = 0, expressed 
in cylindrical coordinates, 

b' r + r _1 6 r + imr~ 1 bg + ikb z = 0, (5) 

shows that it is enough to solve the induction equation for b r and bg only. 

It has been shown 0, @] that for a stationary flow the magnetic field is generated in a resonant 
layer located at r = ro, where the magnetic field lines are aligned with the shear and thus minimize 
their diffusion. In equation (|2|) we see that since there is the same time-dependent factor F(t) 
multiplying angular and axial velocities, this radius is independent of time, and so this surface is 
fixed and given by 

mfi'(ro) + kV'(r ) = 0. (6) 

For more complex time-dependence ro may vary with time, leading to a succession of growing and 
damping magnetic field states 0] • In section IIIII of this paper we will consider the velocity field in 
([2]) with ro fixed, and we assume that ro lies in the fluid, with < ro < 1 (otherwise modes are 
strongly damped). For the more general case of distinct time-dependence for angular and axial 
flows, where ro does vary with time, the equations are set out and discussed in section IIVI 

With a given resonant surface r = ro the leading action of the flow on the field is simply 
advection of the mode by the angular and axial velocities, on the fast advective timescale t = 0(e). 
We take this out of consideration by defining the phase (p(t) as 

<f>(t) = -e- 1 (mO(ro) + kV(r )) [ F(t')dt>, (7) 
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to leave behind only evolution through dynamo action, diffusion and reconnection, on slower time- 
scales. 

Introducing ©, (0} and © in ([I]), we find 

[ed t + (imfi(r) + ikV(r) - imO(r ) - ikV(r )) F(t)]b r 

= e ((L - r" 2 ) b r - 2imr-\) , (8) 

[ed t + {imtt(r) + ifcV(»-) - imO(r ) - ifcV(r )) 

= rQ'{r)F(t)b r + e ((L — r~ 2 ) 6 e + 2imr" 2 6 r ) , (9) 

where the Laplacian operator L is defined by 

L = d? + r- 1 d r -r- 2 m 2 -k 2 . (10) 

In the highly conducting limit e<lwe adopt the smooth Ponomarenko dynamo scaling 0], 

m = e _1 %, & = e~ 1/3 K, r = r + e 1 / 3 ^ t = e 2 / 3 T, (11) 

where r is a time-scale on which the dynamo mode grows, intermediate between the O(l) diffusive 
time-scale and the 0(e) advective time-scale. This scaling gives a magnetic mode localized at the 
radius r = r$ and it is known that the final formulae obtained with this choice of scaling give 
the 'richest' asymptotic picture including both the case m,k = 0(1) and the peak growth rates, 
achieved at m, k = 0(e~ 1 ^ 3 ). Setting 

b r (r,t) =e 1 / 3 b r0 (s,T) + --- , b e {r,t) = &«,(*, r) + •• • , F(t) = F{t), (12) 

together with the 0(r) and V(r) expansion at r = ro, 

n(r) = fi(r„) + e 1/3 ^'(r ) + \e 2 l 3 s 2 tt" \r ) + ..., (13) 
V{r) = V(r ) + e 1 / 3 sV(r ) + ^ 3 s 2 V"(r ) + ..., (14) 

we obtain from (|8|) and ([9]) at leading order in e, 

[0 T + c + ic 2 s 2 F(r) - d 2 ]b r0 = -2iMr - 2 6 eo , (15) 

[d T + c + zc 2 s 2 F(r) - a 2 ]6, = r n'(r )F(T)b r0 , (16) 

with 

c = r ~ 2 M 2 + £T 2 , c 2 = \ (Mfi"(r ) + tfV(r )) . (17) 

From (|15|) and (|16|) we immediately see how the dynamo works: the differential rotation O'(ro) 
stretches radial field b r Q to generate bgo, and the diffusion of bg$ i n cylindrical geometry then 
regenerates 6 r o- For the flows considered below c 2 > 0, and for simplicity we will take this to be 
the case in what follows: there are insignificant changes if this quantity is negative. 

These equations are solved by an exact ansatz involving time-dependent, complex Gaussian 
functions. We put these in, at the same time rescaling using constants a r , ag, a T and a^, to 
eliminate as many parameters as possible, with 

b r o(T,s) \ , \ ( a r f(f) \ n ut^J, 



r, \ I exp(c r) = r ) ' exp(-a h h(T)s ), (18) 
b9o{T,s) J \a 9 g{T) 



where r = a T r is a new time scale. We now fix the constants with 

-2 1/2 „ _ „ _ 1/2 



a r = 2M, ag = r c 2 , a T = a-h = c 2 , (19) 
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and we are left with the following system of ODEs in r to solve 

d f h + 4h 2 = iF(f), 

d f f + 2hf = -ig, (20) 
d f g + 2hg = -$>F(f)f, 

where the constant 'Ss = —2MQ'(ro)/roC2 and the time-dependent factor is F(f) = F(t). (Note 
that modes with more radial structure can be studied, taking the form of two time-dependent 
polynomials times a Gaussian in (|18p . but these will be subdominant). 

The equation for h is nonlinear, while those for / and g are linear. The exponential growth rate 
7 of the magnetic field components / and _g depends only on the parameter 3>, which captures the 
local geometry of the helical streamlines [aj, and the form of the time-dependence F(f). From ([6]) 
and (|17p we have 

9= U n " {ro) (21) 
r \n>(r ) V'(r )J ' 1 ] 

showing that Q depends only on the geometry of the velocity field. The function h(f) gives the 
Gaussian envelope, with Re h > required for exponential localization of the mode. 

So far the system (|20p applies to any time-dependence F(f). Under the rescaling, the two given 
in ([3]) correspond to 

F(f) = cosu)f (ZM), F(f) = 1 + pcoscjf (NZM), (22) 
where the frequencies are linked by 

w = £ - 2/3 c 2 /2 Q = [^e- 1 (mn"(r ) + feU"(r ))] 1/2 Q. (23) 



For the rescaled system (|20p . after a transient, h will become periodic with the same frequency 
u as the forcing F and the linear equations for the magnetic field components (/, g) will take a 
Floquet form: the solution will have exponential growth superposed on periodic behaviour. The 
overall growth rate may be measured as y(@,u) (suppressing p in the non-zero mean case). This 
is linked to the growth rate 7 of magnetic field in the original system © , © (or ([I]) ) with 

7 = e- 2 / 3 [4 /2 7 (^,cu) - c ] = [\e-\wtf\r Q ) + ^"(r ))] 1/2 l(9,G>) ~ r^m 2 - k 2 , (24) 

or, using the definition of £F in (f2T|) . 

7 = [-2mn'(r )/er @} 1/2 7(^,0))- r^m 2 - k 2 . (25) 

In the stationary case F(t) = 1, from ([20]) we find h = ±i 1 / 2 /2 and 7 = -2h ± (i®) 1 / 2 . Then 
from (|24p and taking h with Re h > 0, we obtain the real part of the growth rate as 

Re 7 = e- l/2 {r^ l \m9,\r Q )\) l l 2 - ^e~ 1/2 {\mn"{r ) + kV"{r^)\) l l 2 - rfm 2 - k 2 , (26) 

as previously found 0, 0, [l(| ■ Together with this goes the purely geometrical criterion for Pono- 
marenko dynamo action in highly conducting stationary flow, that |£F(ro)| > 1 at the resonant 
radius r$. Note that formulae (f23"j) - ([26j) . although derived using the scaling (fTTj) . are in fact valid 
for all m and k linked by the resonance condition © . The expansions would give equivalent results 
had we taken m, k = 0(1) and e — > 0, though this would not immediately capture the fastest grow- 
ing modes which are of the scale m, k = 0(e~ 1 ^ s ) as e — > 0. The key assumption in the expansion 
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is that at small e the magnetic field localises in a thin layer. From (jlip and (|18p the width of the 
layer is 

r-r = O (V/V 74 /*- 1 / 2 ) = O (e 1 /4[i( mn //( ro + (ro))] -i/4^-i/2^ (27) 

and this goes to zero as e — > 0; however we should note that this is assuming that Q is fixed as 
we take the limit, so that the magnitude of h is of order unity in the limit. If instead we allow 
u> (or other parameters or wavenumbers) to vary as well, then we need to be careful to check the 
condition that the width given by (|27p is small, to validate the asymptotic theory. For example if 
we fix uj as e — ► we have ui = 0(e 1 ^ 2 uj) — > from (|23|) and as h turns out to be bounded in this 
limit (of low frequencies, so similar to the stationary case) this condition is verified. 

III. RESULTS 

To test the above asymptotic results, we use the flow ([2]) with radial profile Q(r) = 1 — r and 
V(r) = T(l — r 2 ) where T is a helicity factor, for which r$ = —m/2kT and 3> = 4, independent 
of radius. We begin by checking the stationary case, followed by the examples of zero-mean and 
non-zero mean flows in ([3]). 

The growth rate of magnetic field for the asymptotic theory is obtained by simulating f)20f) using 
a fourth order Runge-Kutta scheme, with frequencies and growth rates linked by (f23j) and ([23]) . 
For our given flow we have 

uj = {-kTe- 1 ) l ' 2 Cj, 7 = (-/tre- 1 ) 1 / 2 7-(l + 4r 2 )A; 2 . (28) 

Typically for T = 2 and k = —0.5, we have uj = e~ x I 2 uj and 7 = e~ 1 / 2 7 — 4.25. This growth rate will 
be compared to the one obtained with direct numerical simulation, solving the induction equation 
([1]) without asymptotic approximation, using a Galerkin method for the radial discretization and 
again a Runge-Kutta scheme for the time evolution Note that for the full problem the field 
settles into a Floquet form with B(i + T) = exp(7#T)B(i), where T denotes the period of F(t). 
From dH) and ([7]) it is given by 

7B=1 -i(mn(r )) + kV(r ))T- 1 [ F{t')dt'. (29) 

Jo 

where the phase factor has been reintroduced for a good comparison of frequency measurements. 

A. Stationary flow 

For F(t) = 1 and our given flow, we have 7 = (1 + %)/y/2. From {28} we obtain 73 = 
7 — i(m£l(ro)) + kV(ro)) with 

1 = e~ 1 / 2 ^\k\F (1 + t)- (l + 4r 2 )fc 2 . (30) 

In figure [H the growth rate 75 is plotted against e" 1 for two different values of Vq. The curves show 
a good agreement between both asymptotic and simulation growth rates and frequencies provided 
the magnetic Reynolds number Rm = e _1 is sufficiently large. For e" 1 > 10 3 the difference is 
less than 5% for the growth rate and 0.6% for the frequency. In figure El the modulus of each 
magnetic field component is plotted for several values of e^ 1 . Clearly, increasing e _1 concentrates 
the magnetic field in a thinner layer at ro, and the asymptotic formulation becomes increasingly 




FIG. 1: The magnetic growth rate KejB (left) and frequency ImjB (right) plotted against £ _1 in the 
stationary case F(t) = 1, for m = 1, T = 2 and k = —0.7, ro = 0.35 (black), and k = —0.5, ro = 0.5 (grey). 
The asymptotic solution and the simulation correspond to dashed and full curves respectively. 
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FIG. 2: Modulus of each magnetic field component plotted against r, for (a) e 1 = 500, (b) 1000, (c) 2000, 
(d) 4000, (e) 6000, (f) 10000, for m = 1, T = 2, k = -0.5 and r = 0.5. 

accurate. Denning the layer thickness 5 as the width over which the magnetic energy falls to 
half of its peak value, we find that 5 ~ (9( £ o.27±o.03^ f rom ^ ne simulation. An estimate from the 
asymptotic expressions ([15]) . and using (fTTj) and (fPTj) . leads to 5 = 0{e 1 / 4 ) for m = 1 which is in 
good agreement. However this quantity goes to zero quite slowly with e and so convergence is slow. 



B. Periodic flow with zero mean 

We now consider the second case, with zero mean in the original time-dependence ([3]) or in 
the rescaled version (|22|) . We first work with the asymptotic system and solve (j20|) to obtain the 
growth rate as a function of Q for different values of @, as plotted in figure [3] (left). For 

f < 1 we obtain pure decay Re 7 < 0; for @ = 1.5 the sign of Re 7 depends on cu (positive at small 
u>) whereas for Ql > 2 we find that Re 7 > for all u). Recall that in the stationary case \@\ > 1 is 
necessary and sufficient for dynamo action in the highly conducting limit. 

We can investigate this result further by taking an additional limit of solving the equations (j20[) 
for u) 3> 1. For that we use a new time coordinate u = ujT and set a small parameter Q = u)^ 1 <C 1. 
Then we have, without approximation, from (I20p and dropping the bars to ease notation, 

C l d u h + 4/i 2 = icosu, (31) 
C~ 1 d u f + 2hf = -ig, (32) 
(^dug + Zhg = -Sifcosu, (33) 
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FIG. 3: Zero mean case: growth rate 7 versus ui (left) for, from bottom to top, ^=1,2, 3, 4. The full curves 
give growth rates from integration of ([20]) ; dashed curves and symbols give growth rates from (|40|) . Growth 
rate Re 73 versus (right) for m = 1, V = 2, k = —0.5, r = 0.5 and uj = 200 (black) u = 500 (grey). 
The asymptotic results are shown by dashed curves, simulations by full curves. 

where /, g and h are now functions of u. We set (f(u),g(u)) = exp(fiu)(f*(u), g* (u)) with fi a 
constant Floquet exponent and require /* and g* to be strictly periodic functions of u. We have 

Q^duh + Ah 2 = icosu, 
C 1 »r + C 1 d u f* + 2hf* = -ig\ (34) 
CV<7* + C^duQ* + 2hg* = -^/* cosu. 

Expanding fi, /*, g* and h in powers of £, 

(m, r,g*,h) = (/x, r /») + ccm, r, 5*^)1 + c V r, *o 2 + • • • , (35) 

we solve the system (I34p order by order, using the terms /io, H-i, ■ ■ • to enforce periodicity. This 
leads to 

/* = A + CM + C 2 (^o(2 - 0) cos u + A 2 ) + • • • , (36) 

g* = C@M- sm u ± 2~ 1/2 i) + C, 2 {-@A 1 sinu + £ 2 ) + ••• , (37) 

/i = C(2" 1/2 + isinn) + C 2 C72 + C 3 (C 3 -sin2n + 4\/2icosti) + --- , (38) 

H = C 2 (-2 1 / 2 ± 2" 1 / 2 ^) + • • • , (39) 

where the Ai, Bi and C% are integration constants and we have imposed Re/i > 0. (The Ai are 
arbitrary; the Bi and Cj can be determined in terms of the Ai at higher orders of the expansion.) 
Then the growth rate of / and g is given by 7 = /id) with 

7 ~ uj- 1 2~ 1 / 2 {-2±&). (40) 

This confirms that 7 > only if S> > 2, in the high frequency limit, as seen in figure [3] (left). In 
addition to the results of system (|20p . 7 is plotted versus u) using the asymptotic expansion (|40p . 
We find a good agreement between both, even for moderate values of Q. In figure [3] (right), we 
plot the growth rate -jb for different values of u, from both the asymptotic ODEs (|20p and from 
simulations of the primitive equations (jSJ), ©, showing good agreement provided e~ l is sufficiently 
large. 

Note that formula (|40p indicates a growth rate 7 that increases as the frequency u — ► 0, and 
this perhaps suggests fast dynamo action. The dynamo here would be fast if the growth rate on 



8 



the short, advective time-scale, here given by £7, remains bounded above zero as e — > 0, holding 
the flow fixed. In our flow *2) = 4 which, from (|40l) . leads to 7 ~ \[2jCo and from (f28|) to 

7 ~-v / 2A:re- 1 w- 1 -(l + 4r 2 )A: 2 . (41) 

In the limit of small e, given that —AT > 0, this at first sight appears to be a fast dynamo. This 
formula was derived on the assumption that Co' 1 <C 1, but as e — > for a fixed flow, which includes 
a fixed ui, the assumption becomes violated. As e — > 0, a) — > from (|28p and so we move towards 
the left on figure 3 (left panel): if the asymptotic curves (dashed) continued to grow to the left, the 
dynamo would be fast. However the computed values (solid) saturate for small u and the dynamo 
is in the slow camp, as expected. (As line elements are only stretched linearly the flow fails to have 
Lagrangian chaos, technically positive topological entropy, a requirement for fast dynamo action 
in a smooth flow [la]-) 



C. Critical values for Rm in flows with zero mean 

The asymptotic theory gives an estimate for the critical value of Rm or e for the onset of dynamo 
instabilities, namely from (|24p . 

e x J 2 = Re{[±(mO"(ro) + kV"(r ))} 1/2 l{@,oo)} (rfm 2 + k 2 )~\ (42) 

or from f[25j) . 

e c 1/2 = Re{[-2mO'(ro)/ro^] 1/2 7(^,w)} (rfm 2 + /fc 2 )" 1 . (43) 

For the mode m = 1, the agreement with numerics at onset is poor in figure [T] because the 
critical magnetic Reynolds number is not large enough (contrast the situation in [13]). In the zero- 
mean case the agreement seems to be better as seen in figure [3] (right). In fact we generally expect 
agreement for critical values to improve if there is some other asymptotic parameter to push the 
critical Rm into the small-e, large-Rm regime, provided the condition of thin layer width (|27p is 
satisfied. One possibility could be to take the the limit when the axial flow V(r) is weak or strong 
compared with the angular velocity O(r) (as measured by the helicity factor T). Unfortunately 
direct simulations for r < 1 or T > 1 are difficult to achieve. Instead we consider the limit of 
increasing frequency u. We have already seen that the asymptotic analysis leads to (|4ip provided 
u is large enough. In this case, the threshold should scale as e~ l oc lo and Co tends to infinity 
from (|23p . which is necessary for consistency. Carrying out simulations in order to determine e c 
for values of u in the range 100 < u < 1000, we found that ooe c = 0.209 ± 0.004. This is the 
correct scaling as predicted by (|4ip but the constant is not that predicted, 0.3328. The resolution 
of this paradox is that in this limit with lo = 0(e~ 1 ), we have Co = 0{e' 1 ^ 2 ) — > 00 from (I23p . This 
means that from (|38p h = O(uo~ l ) = 0(e 1 ' 2 ) and so the layer width (|27p is of order unity and does 
not go to zero. The increasing frequency Co is tending to increase the width at the same time as 
the decreasing e is tending to localise the mode, and the two effects cancel out completely. The 
theory gives the correct scaling law but is not asymptotically correct as the magnetic field is not 
localised. This indicates the care that has to be taken with double limits and the usefulness of 
the condition that (|27p be small. (There may be an asymptotic theory appropriate to the limit 
e — ► with elo = 0(1), based on a reduced, finite number of modes in time f similar to those in 
(f36|) ([38]) . but retaining full radial dependence; we leave this for further investigation.) 
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FIG. 4: Non-zero mean case: growth rate 7 versus p (left), for Q = 10. The curves (a), (b), (c), (d) 
correspond to & — 1,2,3,4 respectively. Growth rate Re 75 versus e _1 (right) for m = 1, T = 2, k = —0.5, 
r = 0.5, to = 100 and from bottom to top p = 1,4,8,10. The asymptotic results correspond to dashed 
curves and simulations to full curves. 

D. Periodic flow with non-zero mean 

Here we consider the case with non-zero mean flow, NZM, in ([3]) and (I22p . The growth rate 
7 from the asymptotic ODEs (|20p is plotted against p in figure 2] (left) for u = 10 and different 
values of 0. Taking p = corresponds to the stationary case. Then increasing the fluctuation level 
p may increase the growth rate depending on whether Ql is sufficiently large, the transition being 
for 3l between 2 and 3. This shows that fluctuations may increase the dynamo efficiency (at least 
in the scalings we are using). A different conclusion has been obtained at the dynamo threshold 
which generally increases with the fluctuation rate [3]. When p is increased to large values, the 
mean part of the flow becomes small compared to the fluctuations. We confirmed that this limit 
recovers the results of the previous zero- mean case with appropriate rescaling for 7 and u. 

In figured] (right) the growth rate Re7# is plotted against e" 1 for different values of to. The 
curves show a good agreement between the growth rates from the asymptotic ODEs and from 
simulation of the full system provided e~ l is sufficiently large. The difference is less than 4% for 
e" 1 > 400. 

IV. ANALYSIS FOR TIME- VARYING RESONANT RADIUS 

We now briefly indicate how theory is extended to the more general time dependence, 

U(r,t) = (0,rQ(r)[F(t) + qG(t)],V(r)[F(t) - qG(t)]). (44) 

Here q is a parameter that controls the difference in time-dependence between the axial and az- 
imuthal components and F(t), G(t) are functions of time, of order unity. For example we could 
take a general, single frequency, zero-mean case, 

F(t) = cosujt, G(t) = cos(o>t- (45) 

Now generally the resonant radius becomes a function of time, ro(t), with 

mn\r )(F(t) + qG(t)) + kV'(r )(F(t) - qG(t)) = 0, (46) 

but such variation is found to have a strong damping effect on the field 0] . At the resonant radius 
the shear in the flow is aligned with the helical field lines in the magnetic mode: in the stationary 
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case, as one departs from this radius, the shear changes direction, tending to introduce fine radial 
scales and enhanced diffusion. Moving the resonant radius with a time-dependent flow, then, leads 
to a damping effect on modes, which are strongly suppressed when the field concentration is distant 
from ro (i) (noting that the field cannot readily diffuse in radius to follow ro(t)). 

For these reasons, in our asymptotic framework we will take q to tend to zero in magnitude as 
e — > 0. This makes ro fixed at leading order and we define tq by equation ([6]) as we did originally 
Going through the previous calculations we obtain, in place of () 15f> and (|16p . the equations 

[d T + co + ic lS G(T)+ic 2 s 2 F(T)-d 2 s ]brO = -2iMr-% , (47) 
[d T + c + ic lS G(r) + ic 2 s 2 F(r) - d 2 ]b eo = F(r)r ^ '(r )b r0 , (48) 

with the new, linear term defined by 

G(t) = G(r), q = e 1 / 3 Q, c x = Q(MO'(r ) - KV'(r )). (49) 



Substituting 



M T ' S) ) exp(cor) = ( aJ f \ ) exp(- 0i j(f )* - a h h{f)s 2 ), (50) 
beo\T ,s) J \ a,e9\J) J 



and setting G(t) = G(t), and dj = gives the system, 

d f h + Ah 2 = iF(f), 

d f j + Ahj = iJ2G(f), 

d f f-j 2 f + 2hf = -ig, (51) 

d f -g-?g + 2hg = -9F{f)f. 

We now have a new parameter that quantifies the difference in the time-dependence of azimuthal 
and axial flows, and so the variation in resonant radius, given by 

2 - ci c" 3/4 = ae- 1 ^ m "' (ro) ~ kV ' {ro) (52) 

As for system (I20p . the system (|5ip applies to any time-dependence F(f) and G(f). Again this 
system can be solved numerically to obtain a growth rate. With the specific time-dependence (|45D , 
the growth rate will be a function 7(^,0;,^,$): the phase angle $ quantifies the polarisation of 
the axial and azimuthal components of the flow, in a loose sense. 

In the system (|5ip . we see that changing 2 to — only changes j to —j without affecting 
the other variables. Therefore it is sufficient to consider positive values of 2. Our numerical 
investigations, which we summarize rather than presenting graphically, indicate that compared to 
the curves given in figure [3] for J2 = 0, changing 2 and $ systematically leads to lower values 
of 7, without changing much the shape of the Q and @ dependencies. We can show that 7 is 
7r-periodic in $. We find that 7 is a maximum for $ = 0, a minimum for $ = ir/2 and that 
7($ = vr/4) = 7(3? = 37r/4). Finally 7 is found to be monotonically decreasing with 2 and it 
would be interesting to obtain a proof confirming this observation. 

We can investigate these results further by taking the additional limit of large Q as in section 
EBB. We find that equation ([H holds even for i2 non-zero. In other words the leading order 
growth rate (proportional to w _1 ) is independent of £2 and $ in the limit u — » 00 and is given 
in (|40p . This is very clear in the expansion (|5ip : the j 2 f and j 2 g terms come in at one order 
below what is needed to obtain (|40p . They are asymptotically smaller than hf and hg in the same 
equations, as h and j are both of size £ at leading order. 
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V. CONCLUSIONS 

We extended the theory of the Ponomarenko dynamo in the asymptotic limit of large Rm, to 
the case of a non-stationary flow. We considered only a very simple model for fluctuating flow, but 
one that is nevertheless revealing. Within this class of flows it highlights the basic mechanisms for 
dynamo action, the effect of the motion of the resonant radius in suppressing field generation, and 
the parameter combinations that are relevant at large Rm. Our results include criteria for dynamo 
action at large Rm involving the purely geometrical quantity linked to the rate of change of pitch 
of the velocity shear. For example we find from figure [3] that the geometrical condition \@\ > 2 at 
a given radius is needed for dynamo action with a mode localized there, at high frequencies in the 
zero-mean case for large Rm. 

Note that for stationary flow the corresponding criterion is \&\ > 1 for magnetic field amplifi- 
cation. To see the relevance of these results, consider the spiral Couette flow, which is simply the 
general solution for differentially rotating flow forced by rotating, translating cylindrical boundaries 

fi(r) = Ay + A 2 r~ 2 , V(r) = A 3 + A 4 log r. (53) 

This corresponds to \@\ = 2, and so satisfies the condition in the stationary case: dynamo action 



was observed in [15J]. If the motion of the boundaries is now periodic with zero mean, and sufficiently 
slow that the above functions are just multiplied by F(t) = cos cot, then figure 2 shows that dynamo 
action becomes marginal provided u> is large. Of course the full picture for any boundary forcing 
is complicated by the development of Stokes' layers unless it is slow compared with viscous time- 
scales. Nonetheless the key point is that flows with larger values of \@\ over a range of radii 
are likely to be more efficient as dynamos in non-stationary as well as stationary flows, and this 
consideration could be important in optimizing experiments and understanding experimental or 
numerical results. We also considered flow fluctuations in which components are out of phase, 
leading to a time-dependent resonant radius. The results show that this inhibits the dynamo 



action, confirming previous results obtained for a cellular type of flow [14l ]. 

The fluctuations considered in this paper have a very simple structure, and it would be a 
natural extension to consider fluctuating flows that carried field across streamlines, in other words 
depending also on 9 and z. This would lose the separation of variables employed here and make the 
study more numerical, unless averaging is done analytically, which would generally give an alpha 



effect 16j. Other possible interesting developments would be to consider random time-dependence, 



which increases the complexity of the system, especially close to the threshold 17j, and to include 



some of the effects of nonlinear feedback on the flow field ll[ . 
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